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ABSTRACT 

Transverse oscillations and propagating waves are frequently observed in threads 
of solar prominences/filaments and have been interpreted as kink magnetohydrody- 
namic (MHD) modes. We investigate the spatial damping of propagating kink MHD 
waves in transversely nonuniform and partially ionized prominence threads. Resonant 
absorption and ion-neutral collisions (Cowling's diffusion) are the damping mecha- 
nisms taken into account. The dispersion relation of resonant kink waves in a partially 
ionized magnetic flux tube is numerically solved by considering prominence condi- 
tions. Analytical expressions of the wavelength and damping length as functions of 
the kink mode frequency are obtained in the Thin Tube and Thin Boundary approxi- 
mations. For typically reported periods of thread oscillations, resonant absorption is 
an efficient mechanism for the kink mode spatial damping, while ion-neutral collisions 
have a minor role. Cowling's diffusion dominates both the propagation and damping 
for periods much shorter than those observed. Resonant absorption may explain the 
observed spatial damping of kink waves in prominence threads. The transverse in- 
homogeneity length scale of the threads can be estimated by comparing the observed 
wavelengths and damping lengths with the theoretically predicted values. However, 
the ignorance of the form of the density profile in the transversely nonuniform layer 
introduces inaccuracies in the determination of the inhomogeneity length scale. 

Subject headings: Sun: oscillations - Sun: filaments, prominences - Sun: corona - 
Magnetohydrodynamics (MHD) - Waves 



1. INTRODUCTION 



Waves and oscillatory motions ar e frequently reported in t he observations of solar promi 
nences and filaments (see reviews by lOliver & Ballesterl 12002c iBallester! 120061 : Engvoldl 12008 
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Mackay et al.ll2010|) . In high-resolution observations, the prominence fine structures (threads) of- 
ten display transverse o scillations of small amplitude (e.g., Lin et al]|2005 , 2007 . 2009 : Okamoto et al. 



20071: iNing et a 



(e.g jDiaz et al 



. 2009J ), which have been interprete d as kink magnetohydrodynairnc (M HD) waves 



20021 : bvmova & Rudermanll2005l : iTerradas et al.ll2008l : iLin et al.l l2009). The ob 



served threads in Her images are between 3, 000 km and 28,000 km long, and between 100 km 
and 600 km wide (ILinll2004l : ILin et al.ll2008|) . The threads outline part of much larger magnetic 
flux tubes which are probably rooted in the solar photosphere. The majority of observed periods 
of transverse thread oscillations roughly range between 1 min an d 10 min, but a few detections 
of longer periods of about 20 min have been also informed (e.g., lYi et al J 1 1991c iLinl 12004)) . The 
wavelengths a re usually between 70 km and 8, 000 km, al t hough values up to 250,000 km have 
been reported (lOkamoto et alj|2007|) . Recently, ISoler et al.l (I2010a|) pointed out that the short pe- 
riods and wavelengths are consistent with an interpretation in terms of propagating waves, while 
periods larger than 10 min and wavelengths longer than 100,000 km could correspond to stand- 
ing oscillations of the whole magnetic tube. In the case of standing oscillations, the value of the 
wavelength is not strongly influenced by the thread properties but is mainly determined by the total 
length of the magnetic tube, since the fundamental kink mode wavelength is twice the length of 
the tube, approximately (see details in lSoler et al.ll2010al) . Although there are no direct measure- 
ments of the length of prominence magnetic tubes, this parameter is estimated around 10 5 km. 
This rough estimation is in agreement with the wavelengths reported by lOkamoto et al.l (I2007|) . In 
addition, a common feature of the observations is tha t the oscillations are strongly damped (e.g., 
Terradas et allbooi king et alJlioojt Lin et alJliooi) . 



Motivated by the observational evidence, great effort has been recently devoted to the theoret- 
ical study of both temporal and spatial damping of MHD waves in prominence plasmas. Temporal 
damping is investigated for waves with fixed wavelength, while spatial damping is studied for 
propagating waves with fixed frequency. Both phenomena have been extensively investigated in 
unbo unded and homogeneous prominence plasmas by assuming different damping mechan isms 



Carbonell et al. 



2004, 



20061 boid&orteza et al.ll2007Ll2008h. The reader is referred to blivei 



J2009b . lArregui & Ballested(|2010() . and references therein for a complete account of the theoretical 
works. 

In the case of prominence thread oscillati ons, works so fa r have focused on tempor al damping 



by m echanisms as, e.g., non-adia batic effects (ISoler et al.ll2008l), ion-neutral co ll isions (ISoler et al 



2009bh . and resonant absorption (lArregui et al.ll2Q08L bold : ISoler et alJl2009al B. l201o3) . The con- 
clusions of these works indicate that resonant absorption is efficient enough to provide realistic 
kink mode damping times consistent with the reported strong damping, whereas non-adiabatic 
effects are negligible and ion-neutral collisio ns are only important for shorter wavelengths than 
those observed. In the case of spatial damping. [Pecseli & Engvold (2000) studied the effect of ion- 
neutral collisions but restricted themselves to Alfven waves and kink modes were not investigated. 



-3- 



Although spatial damping of kink wave s has been studied in the context of coronal loops (e.g., 
Pascoe et al.l 12010c iTerradas et al.ll2010ah . to our knowledge no detailed investigation taking into 
account the peculiar properties of prominence threads can be found in the existing literature. The 
purpose of this paper is to fill this gap in the literature, as the recent observations of wave damping 
in solar prominences need to be understood. 

Here, we study the spatial damping of propagating kink waves in prominence threads. Our 
model is composed of a cylindrical magnetic flux tube with partially ionized prominence plasma, 
representing a thread, surrounded by a fully ionized coronal environment. The thread is non- 
uniform in the transverse direction. Resonant absorption and ion-neutral collisions are assumed as 
the damping mechanisms. We use the f3 = approximation, with fi the ratio of the gas pressure 
to the magnetic pressure, and the Thin Boundary approach to describe the effect of resonant ab- 
sorption in the Alfv en continuum using the connection form ulas for the perturbations across the 
resonant layer (e.g. JSakurai et al.lll99ll : lGoossens et al.lll992|) . We determine the dominant damp- 
ing mechanism and obtain analytical expressions for the wavelength, the damping length, and their 
ratio as functions of the kink mode frequency. 

This paper is organized as follows. Section [3] contains the description of the model configu- 
ration and the basic Equations. First, the problem is attacked analytically in Section[3]by adopting 
the thin tube approximation. Later on, the full dispersion relation is numerically solved and a para- 
metric study of the wavelength and damping length of the kink mode as functions of the period is 
performed in Section |4] Finally, the conclusions of this work are given in Section |5J 



2. MODEL AND DISPERSION RELATION 



X 




Y 



► 



Fig. 1 . — Sketch of the prominence thread model adopted in this work. 
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The equilibrium configuration is composed of a straight magnetic cylinder of radius a em- 
bedded in a homogeneous environment representing the coronal medium (see Figured)). We use 
cylindrical coordinates, namely r, ip, and z for the radial, azimuthal, and longitudinal coordinates, 
respectively. The magnetic field is uniform and along the axis of the cylinder, B = Be z , with B 
constant everywhere. Hereafter, subscripts p and c denote prominence and coronal quantities, re- 
spectively. The density within the prominence thread is denoted by p p , while the coronal density 
is p c . Both p p and p c are homogeneous. A transverse transitional layer is included in the radial 
direction, where the density varies continuously between the internal and external densities. We 
do not specify the form of the density profile at this stage. The transverse inhomogeneous length 
scale in the transitional layer is given by the radio l/a, with / the thickness of the layer. This ratio 
ranges from l/a = if no transitional layer is present, to l/a = 2 if the whole tube is radially 
inhomogeneous. Due to the presence of the transverse transitional layer, the kink mode is reso- 
nantly coupled to Alfven continuum modes. The resonance leads to the kink mode damping as the 
energy is transferred to Alfven modes at the Alfven resonance position. This mechanism is known 
as resonant absorption. 

The promin ence plasma is partially ionized and we adopt the single-fluid formalism (e.g., 
Braginskiilll965 ). The ionization degree is arbitrary and is denoted here by the mean atomic weight 
of the prominence material, fx v . This parameter takes values in the range 0.5 < fx v < 1, where 
/ip = 0. 5 corresponds to a fully ionized plasm a and fi v = 1 to a fully neutral gas (see details 
in, e.g., Forteza et al.1 120071 : ISoler et al.ll2009b|) . The external coronal medium is assumed fully 
ionized. The basic MHD Equat i ons governing a pa rtially ionized plasma can be found in, e.g., 
Forteza et al.l (|2007|) : |Pinto et al.l (I2008|) : ISolerl (|2010l) . By assuming the/3 = approximation and 
linear perturbations from the equilibrium state, the basic MHD Equations discussed in this work 
are 

3\ 1 

(1) 



piL = - (V x b) x B ' 

Ot jJ, 



* = V x (v x B) + V x ||| [(V x b) x B] x bJ 



(2) 



where p is the local density, and v = (v r , v v , v z ) and b = (b r , b v , are the velocity and the magnetic 
field perturbations, respectively. Note that v z = in the 8 = approximation. In a partially ionized 
plasma, the induction equation contains a term accounting for Cowling's diffusion, i.e., the second 
term on the right-hand side of Equation Q. Cowling's diffusion represents enhanced magnetic 
diffusion caused by ion-neutral collisions, which is several orders of magnitude m ore efficient tha n 
classical Ohm's diffusion and is the dominant effect in partially ionized plasmas Jcowling 1956 ). 
For this reason, here we neglect Ohm's diffusion and other terms of minor imp ortance present in the 
generalized induction equation (see, e.g., Equation (14) of Forteza et al. Il2007h . Cowling's diffusion 
coefficient, j] c , depends on the ionization degree through p p as well as on the plasma physica l 
conditions. The expression of rj c for a hydrogen plasma can be found in, e.g.. |Pinto et al.1 (I2008|) 
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and Sol er et al.1 (|2009b() . whereas for a plasma composed of hydrogen and helium see lSoler et al.. 
(i2010bh . As the effect of helium is negligible for realistic helium abundances in prominences, here 
we consider a pure hydrogen plasma. The effect of Cowling's diffusion is neglected in the external 
medium because the corona is assumed fully ionized. 

We follow an approach based on normal modes. Since (p and z are an ignorable coordinates, 
the perturbations are expressed proportional to exp (imcp + ik z z - icot), where co is the oscillatory 
frequency, k z is the longitudinal wavenumber, and m is the azimuthal wavenumber (m = 1 for the 
kink mode). Alternatively, the problem could b e investigated by means of time-dependent simula- 
tions of driven waves as in Pascoe et al. However, in the linear regime the different values 
of m and k z are decoupled from each other, and a normal mode analysis is a simpler procedure for 
linear waves. If Cow ling's diffusion is neglected, our configuration corresponds to that studied by 
Terradas et al] (I2010a|) for propagating kink waves in coronal loops. We extend their investigation 
by incorporating the effect of Cowling's diffusion due to ion-neutral collisions 

By using the Thin Boundary (TB) approach (see details in, e.g., Goossens et al.l2006 : Goossens 



20081) . the analytical dispersion relation for resonantly damped kink waves prop agating in a trans 
versely nonuniform and partially ionized prominence thread was obtained by ISoler et al.l (12009c . 
Equation (25)). If partial ioniza tion is not consid ered and the effect of Cowling's diffusion is ab 



sent, th e dispersion relation of ISoler et al] (|2009ch redu ces to that investiga ted by iTerradas et al 



(I2010al Equation (28)) for kink waves in coronal loops. ISoler et al] (I2009q) checked that the so- 
lutions of their dispersion relation are in excellent agreement with the solutions obtained from the 
full numerical integration of t he MHD equations beyond the TB approximation. Therefore, the 
dispersion relation derived by ISoler et al] (|2009d) correctly describes the kink mode behavior in 
ou r model and compli cated numerical integrations are not needed. The dispersion relation obtained 
by lSoler et al] (l2QQ9q) in the case of a straight and homogeneous magnetic field is 



n c 



K' (n c a) 



J'm ( m P a ) 



= -17T- 



m 2 /r 2 A 
' oj 2 \d r p\ r 



m 



Pc (co 2 - k 2 z v\ c ) K m (n c a) Pp _ j m 

whe re J m and K m are the Bessel function and the modified Bessel function of the first kind of order 
( Abramowitz & Stegun 1972|) . respectively, and the quantities m p and n c are defined as 

(c 2 -k 2 Y 2 Av ) (ky Ac -co 2 ) 



(3) 



r 2 

1 Ap 



(4) 



'Ac 



where T 2 Ap 
with v Ap = 



'Ap 



- icorjc is the modified prominence Alfven speed squared (|Forteza et al]|2008|) . 
and v Ac 



^ -j= the prominence and coronal Alfven speeds, respectively, and 

// = An x 10~ 7 N A~ 2 the magnetic permeability. In addition, r A is the Alfven resonance position, 
and \d r p\ rA is the radial derivative of the transverse density profile at the Alfven resonance position. 
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Soler et al.1 (I2009q) studied the temporal damping of kink waves, hence they assumed a fixed, 



real k z and solved Equation © to obtain the complex frequency. Here, we investigate the spatial 
damping and proceed the other way round, i.e., we fix a real co and solve Equation © to obtain the 
complex wavenumber. Then, the period, P, wavelength, A, and damping length, L D , are computed 
as follows, 

2tt 2n 1 

P = — ^=7—' L D = —, (5) 

CO KrR K z l 

with fc zR and k A the real and imaginary parts of k z , respectively. 



3. ANALYTICAL APPROXIMATIONS 

Some analytical progress can be performed before solving Equation © by means of numerical 
methods. To do so, we adopt the Thin Tube (TT) limit, i.e., A/a » 1. A fist-order expansion of 
Equation © gives the dispersion relation in both the TT and TB approximations, namely 

I 2 V 2 T 2 \ , I 2 , 2 2 \ ■ m PpPc (^ 2 -^ r Ap)(^-^Ac) ... 

Pp(^ -^r Ap )+p c (o, -k z v Ac )-m-j^- =0. (6) 

If both Cowling's diffusion and resonant absorption are omitted, the solution to Equation © is 

kl = - r , (7) 

with c 2 = \ the kink speed squared. Equation © corresponds to the ideal, undamped kink 

mode. The solutions to Equation © considering the different damping mechanisms are discussed 
next. 



3.1. Damping by Cowling's diffusion 

In the absence of transverse transitional layer, i.e., l/a = 0, resonant absorption does not take 
place and the damping is due to Cowling's diffusion exclusively. In such a case, the third term on 
the left-hand side of Equation © is absent. We write the wavenumber as k z = k zR + ik zi and use 
Equation © to obtain the exact expressions for k 2 R and k^, namely 



k 2 



*5 



2„2 



CO c 



4 , 7-2 



1 

2 ct + co 2 fji 



co z ni 

'1 + — £ + 1 



to 2 rrL 
+ — £-1 



(8) 



(9) 
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with r/c = -^—J]c. By combining Equations ([8]) and ©, we compute the ratio of the damping 

Pp+Pc 

length to the wavelength as 

L» = k^ = 1 c l + S + 
A 2nk a 2n cofjc 

Equations (l8l)- ([T0l) are exact expressions that can be further simplified depending on the value 

2-2 2-2 

of the ratio For — ^ «; 1, i.e., in the limit of low frequency {to small) and/or large ionization 
degree (77c small), Equations (ISl)-(fTUl) simplify to 



-.2 / 1 _l ^ \ C£ 



1 4-2 1 4-2 

2 l CO T c l OTT^ 

" ?;^7^p— (l2) 

l ^ 



/I n corjc 



(13) 



On the contrary, if ^-fi » 1, i.e., high frequency and/or small ionization degree, the equivalent 



c, 



expressions are 



T n 1 tO 

Zc k / 7c 

i * In 



(15) 



.2=2 



Thus, for ^-^ ^ 1 the ratio of the damping length to the wavelength is inversely proportional to 
both to and fjc, and k^ R coincides with the ideal value (Equation ©). This case corresponds to a 

2-2 

weakly damped kink mode. On the other hand, for » 1, Ld/A is independent of to and i) c and 
the wave behavior is governed by diffusion. By assuming typical values for the parameters in the 
context of oscillating prominence threads, e.g., P = 3 min, B = 5 G, and p p /p c = 200, we obtain 

2-2 2-2 2-2 

—3^ «6x 10~ 17 for /i p = 0.5 and — ^ « 1.6 x 10~ 4 for /2 p = 0.99, meaning that the case — ^ «c 1 

k k k 

is more realistic in the context of oscillating threads even for an almost neutral plasma. 



3.2. Damping by resonant absorption and Cowling's diffusion 



Next, we take the case l/a ± into account and study the combined effect of resonant absorp- 
tion and Cowling's diffusion. The third term on the left-hand side of Equation © is now present. 
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As before, we write k z = k Z R + ik z \ and put this expression in Equation ©. Since it is very difficult 
to give exact expressions for k zR and k d in the general case, we focus on L D /A and restrict ourselves 

2 -2 I 1 I 1 

to ^ 1- Following the procedure of iTerradas et al.1 (I2010ah . we assume weak damping, i.e., 

C k 

k z i «; k zR , and negle ct terms with k? v Th e following process is long but straightforward, and we 
refer the reader to ITerradas et al.1 (|2010ah for details. Finally, we arrive at the expression for the 
ratio of the damping length to the wavelength as 



Ld 
A 



U)T]c 



n- 



m I Pp-Pc 
+ 



T ap p +p c 



(16) 



where the first term within the parentheses accounts for Cowling's diffusion and the second term 
for resonant absorption. The factor T in the second term takes different values depending on 
the density profile wi thin the inhomogeneous layer. For example, T = 4/ n 2 for a linear profile 
(Go ossens et al.ll2002|) . while T = 2/n for a sinusoidal profile with r A ~ a (|Ruderman & Roberts 



20021) . If the term related to resonant absorption is absent, Equation (fT6l) reverts to Equation (fT3l) . 
On the other hand, if the term related to Co wling's diffusion is dropped, Equation (TT6l) coincides 
with Equation (13) of ITerradas et al.1 (I2010ah . 



The relative importance of the two terms in Equation (1161) can be assessed by performing their 
ratio as 

(LdM)rA . _n- a UT]C Pp + Pc 0)T]c Pp 



6 = 



l elm p p - p c 



/ elm Pp - p c 



(17) 



where (Lu/A)^ and (L D /A) C stand for the damping ratio by resonant absorption and Cowling's 
diffusion, respectively. By considering as before P = 3 min, B = 5 G, and p p /p c = 200, and 
adopting a linear profile with l/a = 0.2, we obtain e « 8 x 10~ 8 for /i p = 0.5 and e « 0.12 
for jlp = 0.99, meaning that in the TT limit resonant absorption dominates the kink mode spatial 
damping for typica l parameters of thread oscillations. This result is equivalent to that obtained by 
Soler et al.l (|2009q) in the case of temporal damping. 



4. NUMERICAL RESULTS 

Now, we solve the dispersion relation (Equation ©) by means of standard numerical proce- 
dures. In the following figures, both the wavelength, A, and the damping length, L D , are plotted 
in dimensionless form with respect to the thread mean radius, a. The period, P, is computed in 
units of the internal Alfven travel time, r Ap = «/va p - Unless otherwise stated, the results have 
been computed with p p = 5 x 10~ n kg m 3 , p p /p c = 200, and B = 5 G. With these parameters, 
v Ap ~ 63 km s" 1 and, for a = 100 km, r Ap « 1.59 s. 

Figure [2] displays A/ a, L D /a, and L^/A versus P/r Ap for the case l/a = 0, i.e., the damping 
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is due to Cowling's diffusion exclusively. We compute the results for different values of p p . The 
shaded areas in Figure [2] and in the other figures represent the range of observed periods of trans- 
verse thread oscillations, i.e., 1 min - 10 min, corresponding to 40 < P/t Av < 400, approximately. 
Regarding the wavelength, we see that the effect of Cowling's diffusion is only relevant for periods 
much shorted than those observed. This is in agreement with Equations (fTTT) and (fT4"l) . On the other 
hand, an almost neutral plasma, i.e., p. v — » 1, has to be considered to obtain an efficient damping 
and to achieve small values of L D /A within the relevant range of periods. Although we do not 
know the exact ioniz ation degree in prominence threa ds, such very large values of fi v are probably 



unrealistic (see, e.g., iGouttebroze & Labrossell2009i) . The analytical expressions for A, L D , and 



L D //l in the TT case given by Equations ©, ©, and (flOl) . respectively, are in good agreement with 
the full results in the whole range of periods (see symbols in Figure [2]). 

Next, we study the general c ase l/a # 0. We adopt a sinusoidal density profile within the 



inhomogeneous transitional layer (|Ruderman & Robertsll2002l) . As the Alfven resonance position 



r A , is needed for the computations of the resonant damping, we follow a two-step procedure. First, 
we solve the dispersion relation for a fixed to in the case l/a = and determine k Z R. Then, we 
assume that the value of k zR is approximately the same in the case l/a + 0, meaning that the 
resonant condition is co = k z ^y A (t"a)- In the case of a sinusoidal profile, the expression of the 
resonant position can be analytically obtained from the resonant condition as 

' Pp + Pc 2v Ap k zR Pp 



I 

r A = a + - arcsin 
n 



{p p - p c co 2 p p - p c 



(18) 



Finally, we compute \d r p\ rA using the previously determined r A by means of Equation (fT8l) and 
solve the dispersion relation with these parameters to obtain the actual k zR and k z \. Figure [3] shows 
the results of these computations for different values of l/a when the ionization degree has been 
fixed to jlp = 0.8, whereas Figure H]displays the equivalent computations for different values of p p 
when the transverse inhomogeneity length scale has been fixed to l/a = 0.2. Since the wavelength 
is not affected by the value of l/a and has the same behavior as in Figure (2a), both Figures [3] and 
|4] focus on L D /a and L D /A. We obtain two different behaviors of the solutions depending on the 
period. For small P/r Ap , the damping length is independent of l/a and is governed by the value 
of pp. On the contrary, for large P/r Ap the damping length depends on l/a but is independent of 
// p . This result indicates that resonant absorption dominates the damping for large P/r Ap , whereas 
Cowling's diffusion is more relevant for small P/r Ap . The approximate transitional period, namely 
P tr , in which the damping length by Cowling's diffusion becomes smaller than that due to resonant 
absorption can be estimated by setting e ~ 1 in Equation (fTTT) and writing P tr = In /to. Then, one 
obtains 

P b ~ In T 1 —, = 2n T-r— 7 . (19) 

/ c^mpp - p c / cpnp p - p c 

This transitional period is in good agreement with the numerical results (see the vertical dotted line 
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Fig. 2. — Results for the kink mode spatial damping in the case l/a = 0: (a) A/a, (b) L D /a, and (c) 
Ld/A versus P/r Ap for ju p = 0.5, 0.6, 0.8, and 0.95. Symbols in panels (a), (b), and (c) correspond to 
the analytical solution in the TT approximation given by Equations ([8]), ©, and (flOl) . respectively, 
while the horizontal dotted line in panel (c) corresponds to the limit of Ld/A for high frequencies 
(Equation (fTBT)). The shaded area denotes the range of observed periods of thread oscillations. 
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Fig. 3. — Results for the kink mode spatial damping in the case l/a ^ 0: (a) L D /a and (b) L D /A 
versus P/t Ap for l/a = 0.05, 0.1, 0.2, and 0.4, with /} p = 0.8. Symbols in panel (b) correspond to 
the analytical solution in the TT approximation given by Equation (fT6l) . while the vertical dotted 
line is the approximate transitional period given by Equation ([T9l for l/a = 0.1. The shaded area 
denotes the range of observed periods of thread oscillations. 




Fig. 4. — Results for the kink mode spatial damping in the case l/a ± 0: (a) L D /a and (b) L D /A 
versus P/r Ap for ju p = 0.5, 0.6, 0.8, and 0.95, with l/a = 0.2. Symbols in panel (b) correspond to 
the analytical solution in the TT approximation given by Equation (fT6l) . The shaded area denotes 
the range of observed periods of thread oscillations. 
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in Figure [3tb)). In addition, we see that P tr is much smaller than the typically observed periods, 
indicating that resonant absorption is the dominant damping mechanism in the relevant range. 

Finally, we check that the analytical approximation of L D //l given by Equation (fT6b provides 
an accurate description of the kink mode spatial damping in the relevant range of periods (compare 
the symbols and the solid lines in Figures [3jb) andHJb)). 



5. DISCUSSION AND CONCLUSION 



In this paper, we have studied the spatial damping of kink waves in prominence threads. 
Resonant absorption and Cowling's diffusion are the damping mechanisms taken into account. 
Both analytical expressions and numerical results indicate that, in the range of typically observed 
periods of prominence thread oscillations, the effect of Cowling's diffusion (and so the ionization 
degree) is negligible. On th e other hand, r e sonant absorption provides an efficient damping in 
agreement with the study of iTerradas et al.1 (I2010a|) i n the context of co ronal loop oscillations. 
These conclusions are equivalent to those obtained by ISoler et al.1 (I2009ch in the case of temporal 
damping. 

We point out that small values of L D /A are obtained by resonant absorption in the observa- 
tionally relevant range of periods, which is consistent with the reported strong damping of the 
oscillations. The analytical estimation of L D /A given by Equation (|T6b is very accurate in the ob- 
servationally relevant range of periods, and the contribution of Cowling's diffusion can be dropped 
from Equation (fT6l) because the plasma ionization degree turns out to be irrelevant for the damping. 
Therefore, for kink modes (m = 1) the radio L D /A simplifies to 



A I p p - p c ' 



which coincides with the expression provided by ITerradas et al.l (|2010a|) . As 



Pp+Pc 
Pp-Pc 



(20) 



1 for typical 



prominence and coronal densities, this factor can be dropped from Equation dzOj), meaning that the 
ratio L D /A depends almost exclusively on the transverse inhomogeneity length scale, I /a, and the 
form of the density profile through T as 



^2 ~ <F- 
A ~ V 



(21) 



In the case of coronal loop oscillations studied by ITerradas et al. ( 2010a ). the factor cannot be 
dropped from their expressions, meaning that in coronal loops the ratio L D /A significantly depends 
on the density contrast. Therefore, information about the parameters l/a and T in prominence 
threads could be determined by using Equation (1211) along with accurate measurements of the 
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damping length and the wavelength provided from the observations. However, since the precise 
form of the transverse density profile in prominence threads is unknown, we have to assume an 
ad hoc profile, i.e., a value of T , to infer the transverse inhomogeneity length scale from the 
observations, which can introduce some uncertainties in the estimation of l/a. 

For example, let us assume that the ratio Ld/A has been determined from an observation of 
damped kink waves in a prominence thread and we want to compute the transverse inhomogeneity 
length scale of the thread. For simplicity, we consider that the transverse density profile in the 
inhomogeneous layer is either linear or sinusoidal. Denoting as (l/a\ in the value of l/a computed 
assuming a linear profile, and (l/a) sin the corresponding value for a sinusoidal profile, the relation 
between both of them is 

^ = = - 1-57, (22) 

pointing out that the relative uncertainty of l/a is larger than 50%, and the inaccuracy could be even 
larger if other profiles are considered. This fact should be taken into account in future seismological 
determinations of this parameter. 

The present investigation is a first step for the study of the spatial damping of kink waves 
in prominence fine structures. Here, we have adopted a simple model of a prominence thread. 
Some effects that might influence the kink mode propagation and damping are not included in 
the present paper. Among them, plasma inhomogeneity along the thread may affect somehow 
the amplitude of a propa gating kink mode, wh ereas the presence of flows affects the damping by 



resonant absorption (see lTerradas et al.ll2010b|) . The influence of these and other effects will be the 



subject of forthcoming works. 
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